Directional 'superradiant' collisions: bosonic amplification of atom pairs emitted from 

an elongated Bose-Einstein condensate 
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We study spontaneous directionality in the bosonic amplification of atom pairs emitted from 
an elongated Bose-Einstein condensate (BEC), an effect analogous to 'superradiant' emission of 
atom-photon pairs. Using a simplified model, we make analytic predictions regarding directional 
effects for both atom-atom and atom-photon emission. These are confirmed by numerical mean- 
field simulations, demonstrating the the feasibility of nearly perfect directional emission along the 
condensate axis. The dependence of the emission angle on the pump strength for atom-atom pairs 
is significantly different than for atom-photon pairs. 

PACS numbers: 03.75.Fi, 03. 75. Be, 03.75.-b 



One of the most intriguing nonlinear/quantum phe- 
nomenon observed in an atomic BEC to date has been 
the observation of spontaneous directionality in Rayleigh 
scattering of laser hght from an elongated condensate 
101 . In this experiment a pair of condensate momentum 
sidemodes were generated along with a pair of superra- 
diant light pulses in a cigar-shaped condensate, whereas 
an isotropic burst of non-condensate atoms and light was 
observed from a spherical BEC. 

From a quantum-optics point of view, the laser-driven 
BEC can be viewed as a coherent source of correlated 
atom-photon pairs [||. Various schemes have been pro- 
posed recently in which a BEC emits correlated atom- 
atom rather than atom-photon pairs. One such scheme 
makes use of microwave stark shifts to manipulate 
the kinematics of spin-mixing collisions in a spinor BEC. 
Other schemes involve the dissociation of a BEC of di- 
atomic molecules . Provided a pure diatomic molecu- 
lar condensate is eventually realized P-p^, its dissocia- 
tion into atom pairs could be used to generate various 
macroscopic entangled quantum states A third, 

currently feasible, scheme proposes generating correlated 
atom pairs by stimulated Raman coupling via a molecu- 
lar state . 

From both fundamental and applied considerations it 
is important to know whether or not spontaneous di- 
rectionality is achievable in the case of atom-atom pair 
emission from an elongated BEC. If a single pair of new 
condensates could be generated in this manner, than by 
virtue of the atoms having been emitted as correlated 
pairs there would be strong non-classical correlations, i.e. 
'squeezing', between the new condensates, with impor- 
tant potential applications such as sub-shot-noise atom 
interferometry p^ ]. 

While directional effects were postulated in ||^, the 
analysis there relies on a Markovian approximation in 
which a portion of the atomic field inside the BEC is 
treated as a reservoir with an infinitely short 'memory'. 
This approach, while justified for a light field inside the 



BEC g], can not be applied to slowly-evolving atomic 
fields, whose long 'memory' is responsible for directional- 
ity. Due to this flaw, no conclusions regarding the physics 
of directionality can yet be drawn. We note, however, 
that there is no inconsistency with treating the atomic 
fleld outside the BEC as a Markovian reservoir, because 
once atoms exit the BEC they are irreversibly lost and 
can no longer influence the ensuing physics inside the 
BEC volume. 

The purpose of this Letter is to analyze the physics 
of spontaneous directionality when correlated pairs of 
atoms are generated from an elongated condensate; 
whether via spin-mixing collisions, photodissociation of 
molecules, or molecular-resonance Raman spectroscopy. 
We show that by introducing a large mass ratio of the 
paired atoms we can reproduce the physics of superra- 
diant Rayleigh scattering. In the limit of equal masses, 
on the other hand, we smoothly map onto a new regime 
in which directionality still exists, but with significant 
differences. While generating number-squeezed states is 
an important application of 'superradiant collisions', we 
will not analyze in detail the quantum statistics at the 
present time. 

We consider, for concreteness, a photodissociation 
scheme using two lasers to set up a bound-free Raman 
transition via a far-detuned electronically excited molec- 
ular state. The two-photon resonance is tuned into the 
dissociation continuum, so that excess energy goes into 
the relative motion of the atom pair. Using copropagat- 
ing lasers, it is possible to transfer sufficient energy to 
the relative motion with negligible transfer of center-of- 
mass momentum. As a variation on this scheme we note 
that the initial state can also be an unbound atomic BEC 
state, as in [Q. The coupling via the molecular excited 
state will transfer energy to the relative motion of atoms, 
again achieving the emission of atom-atom pairs from the 
BEC. As experiments to implement are underway us- 
ing a spherical BEC with counter-propagating lasers at 
zero detuning, only minor changes in the laser parame- 
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ters and BEC aspect ratio would be required to observe 
the effects we describe. We note that our model will also 
describe collision-based schemes [^,| 14 1 with only incon- 
sequential adjustments. 

We begin our analysis with the Hamiltonian for atom- 
pair generation: 
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Here, the operator Hi is the free Hamiltonian of the 
quantum field 4'i{v), where i describes different atomic 
states or species, N is the particle number, Q.ij is the 
two-photon Rabi frequency, K is the two-photon re- 
coil momentum, and A, which supplies kinetic energy 
to the relative motion of the atom pairs, is the two- 
photon frequency difference minus the frequency corre- 
sponding to the chemical potential of the molecular state 
and the molecular binding energy. The function <&(r) 
is the molecular BEC wavefunction, normalized to unity. 
For spin-exchange collisions or molecular resonance spec- 
troscopy $(r) would be the square of the atomic conden- 
sate wavefunction. 

The Heisenberg equations of motion for the atomic 
field operators in the rotating frame ((/'i(r) 
0i(r) exp[— i((5/2)i]) read 



m,$(r)e*'^-"-0j(r). 



(2) 



where length and time are rescaled as r ^ Lr and t 
with T = miL^/h, k = LK. and 6 — tA. 

The free Hamiltonian Hi contains kinetic energy, trap 
potential and condensate mean-field contributions. Atom 
pairs generated within the BEC volume could therefore 
be described initially by a superposition of Eigenmodes 
of Hi. The finite energy width of these superpositions 
is then consistent with the finite time scale on which 
the wavepacket freely evolves out of the BEC volume. 
Consequently, the scattering states can be represented 
by a restricted basis-set of eigenstates of a 'virtual cav- 
ity' within the volume of the BEC The losses due to 
translational wavepacket motion can then be modeled by 
a linear loss term. 

Applying this approach, we use eigenstates of a box of 
length L and width L/a (with periodic boundary con- 
ditions as the lack of a physical box precludes stand- 
ing wave states) to expand (piir) — X)q "^^iq, where 
the states are normalized in dimensionless units and the 
summation is over a grid of q values satisfying q = 
a{ix + my) + nz, for any integers m, and n. To treat 
the free atomic Hamiltonians, we then approximate 

- l^-ffi, e"''''CiqJ = 2 '2*'*''" [^'^»q + »7»q] Cjq + /jq (3) 



where Wiq = ^i|qp is the kinetic energy term (/i^ = 
mi/mi). The decay of population from mode q, 
due to translational motion of the atomic wavepacket, 
is generally non-exponential. In order to approxi- 
mate these losses as exponential decay we introduce 
a loss rate proportional to the atomic velocity q 
over the BEC dimension along q according to 7iq = 

1/2 

I3^i [(aq • x)^ + (aq • y)^ -I- (q • z)^] , where the fitting 
parameter (3 ~ 0[1] then depends on the exact shape 
of the BEC wavefunction. Lastly, /^q is a standard 
quantum-Langevin noise operator which must be in- 
cluded in any damped quantum system. 

Approximating the normalized BEC ground state by 
the lowest energy box-eigenstate, we arrive at a finite set 
of coupled Heisenberg equations 
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where Xij = T^ij \fri, n being the condensate density. We 
note that we have dropped the noise operator as we are 
interested in the presence of instabilities rather than in 
their detailed quantum statistics. Due to orthogonality, 
the modes {c^q} are only coupled to the modes {cj(k_q)}. 
Equation (^ is thus a finite set of linearly coupled oper- 
ator equations, whose eigenfrequencies characterize the 
dynamics. The existence of an eigenfrequency with pos- 
itive imaginary part would indicate the onset of para- 
metric instability, in which initial spontaneously emit- 
ted atom pairs are amplified, leading to an exponential 
buildup of a macroscopic mode population. 

In order to demonstrate the connection between atom- 
pair generation and superradiant Rayleigh scattering, we 
consider first the case of a two-species atom field where 
/i2 ^ 1 and Xij = X^ii^j2, corresponding to the disso- 
ciation of a heteronuclear molecule into one heavy and 
one light atom. Under these conditions, the more mas- 
sive atomic field can build up a macroscopic population 
while the less-massive field, owing to its faster escape ve- 
locity, will remain only weakly populated. This situation 
is analogous to the regime of BEC superradiance, with 
the light atoms corresponding to scattered photons. The 
rapid damping of the light-atom field allows for its adia- 
batic elimination, which leads to a single linear equation 
for the operator ciq. Solving this equation, we find that 
the mode occupation (cjqCiq) grows exponentially at the 
rate 
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This expression reproduces the gain expression for super- 
radiant Rayleigh scattering if we replace the atomic 
dispersion and decay terms with those of photons. The 
heavy atoms, therefore, will be ejected at an angle q to 
the long axis satisfying q = k± y/5/ii2Z, as was observed 
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in [Q. Equation (||) shows that there is no saturation of 
directionality with increasing pump strength x- While 
this is essentially the gain expression given in , the va- 
lidity condition, /i2 3> 1, is not satisfied in that system. 

It is interesting to note that while it is the build up 
of macroscopic populations in the slow-moving atomic 
modes which leads to exponential growth and direction- 
ality, it is the losses of the fast-moving atoms or pho- 
tons which is minimized. This state of affairs can be un- 
derstood if we view the light-atom (or photon) 'virtual 
cavity' modes as intermediate states through which the 
particles pass before decaying irreversibly into a reser- 
voir of modes. It is well known in quantum optics that 
the rate of coupling to a continuum is reduced if a short- 
lived intermediate state is introduced, as its energy width 
acts as a detuning making the transition increasingly off- 
resonant as the lifetime decreases. 

In the remainder of this Letter we focus on the the 
photodissociation of homonuclear molecules. Using co- 
propagating lasers, the momentum transfer associated 
with the bound-free transition is negligible compared to 
the momentum spread of the BEC, i.e. k w 0. The 
'virtual cavity' modes then decouple into ciq and ci_q 
pairs, resulting in a quadratic eigenvalue equation, with 
solutions 
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The condition for positive exponential gain (Imjajq} > 0) 
is therefore 
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leading to exponential growth of the mode occupation at 
the rate 
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Equation (0) shows that bose-stimulation leads to ex- 
ponential amplification provided that the pumping is 
sufficiently strong to overcome the combined effects of 
losses and de-phasing due to energy mismatch. Energy 
is conserved on a sphere in momentum space centered 
at q = with radius ^fS. The gain along this sphere, 
given by Gq = 2|x| — 7q, attains a maximum when the 
directionally-dependent loss 7q is minimized. As a result, 
modes directed along the long-axis of the BEC will grow 
the fastest. 

The minimum loss rate value of (3^/5 along the conden- 
sate axis, sets a power threshold of 2x > P\f6 for end-fire 
amplification. In comparison, the stimulation threshold 
for side-fire modes, 2% > ^cx\f6 is larger by a factor of 
the aspect ratio a. We therefore identify three regimes 
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FIG. 1. The exponential instability threshold angle the 
gain feature HWHM angle 9hw, and the single-mode diffrac- 
tion angle 8d as functions of the dimensionless control param- 
eter y. 

as the parameter y — 2x/ iP\/S) is varied with respect to 
these two thresholds : 

Case I: y < 1. When the pump intensity is below 
both axial and radial thresholds, there is simply isotropic 
spontaneous emission of atom pairs. 

Case II: 1 < y < a. Between the thresholds, the 
exponential gain condition is satisfied for modes with an 
emission angle 9e from the z-axis. In general, a random 
spot pattern will appear within O^., with characteristic 
spot-size determined by the angle subtended in k-space 
by a single resonant end-fire mode 9d ~ 4a/ v^. 

Just above the longitudinal threshold {1 < y a) 
we have 9e ~ 1/a. If 6e is smaller than the spot-size 
Od, only one end-fire mode in each direction along the 
condensate axis will grow to macroscopic size. This will 
be the ideal situation for generating a number-squeezed 
pair of condensates. 

For y extremely close to the lower threshold the expo- 
nential growth may be slow enough that losses are not 
negligible during the build-up of the number-squeezed 
state. In this case jets, i.e. out-coupled atom-laser 
beams, will be observed exiting the BEC volume along 
the z-axis. This can lead to a decay of the number- 
squeezing, as the loss of atoms from the end-fire modes is 
an uncorrelated random process. By slightly increasing 
the pump strength slightly these can be made negligible 
without sacrificing the narrow emission angle. 

As the pump intensity increases further into the in- 
termediate regime, both the width of the resonance shell 
and the solid angle of emission increase dramatically with 
y due to power broadening . From Eq. (||), we find 
that the angle with respect to the BEC axis within which 
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amplification can occur is 9^ = 
to mode competition, however, 
place within the angle 9hm = s 
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most emission will take 

n-l [j-^iy + iy2-A 

corresponding to the half-maximum of the gain feature 
of Eq. (||). In Figure 1 we plot 9e and 9hm versus y, 
which shows both saturating at 7r/2 by y = 20. The pri- 
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mary result of this Letter, therefore, is that two-mode 
directionality occurs only in a narrow operating regime 
1 < y < j/cr, where ycr is the critical value where the 
emission angle equals the diffraction angle Of, = 9d In 
contrast, no such narrow limitation exists in directional 
superradiant Rayleigh scattering . 

Case III: y > a: When the pump intensity is above 
both thresholds, the emission angle is comparable to the 
saturation angle 7r/2. While a slight gain differential may 
still give some preference towards smaller angles, in the 
limit y ^ a the gain differential becomes negligible, re- 
sulting again in isotropic emission. 

In order to test the above predictions, based on the 
analysis of Eq. (||), we solve Eq. (||) numerically 
for the dynamics of the atomic field-operator cf). The 
dominant quantum-field effect is amplification of sponta- 
neously emitted atom pairs in the vicinity of the insta- 
bility. Therefore, Eq. (H) may be solved to a very good 
approximation, by propagating its classical (mean-field) 
limit (pi ^ (pi, with an initial atomic field seeded by ran- 
dom noise of order 1/\/N. We note that this approach 
has been widely used to simulate parametric superfluo- 
rescence [ p^ in quantum optics. The increase in emission 
angle with y is clearly seen in the 2-dimensional numeri- 
cal simulations. 

In Figure 2 we plot the atomic density in k-space for 
various values of the control parameter y. The simula- 
tions use a Gaussian condensate wavefunction with as- 
pect ratio a = 10, a detuning (5=10^ and are integrated 
via a split-step FFT algorithm until such time as 10^ 
atoms are emitted. The free Hamiltonian is taken to 
contain only the standard kinetic energy term for sim- 
plicity. Adding a trap potential will have no significant 
effect on the timescales we consider. In Figure 2a we 
show the case y — 2, in which a single pair of end-fire 
modes is observed. The case y = 10 is shown in Fig 
2b, where we see that multiple modes are generated, all 
falling within the HWHM of Eq. (^, as predicted by 
our analysis. In figures 2c and 2d we plot the cases 
y — 2Q and y — 100, respectively, demonstrating the 
saturation of directionality as y is increased beyond the 
radial threshold. The contours in the figures correspond 
to our analytical predictions as described in the caption. 
Quantitatively, the growth rate of the end-fire modes, 
predicted by Eq. to be G = 2x — PVS, was numeri- 
cally found to be G = 1.8% — 5.4-\A, which corresponds 
to /3 « 6 for a Gaussian wavepacket. 

Interestingly, these simulations reveal the appearance 
of correlated pairs of modes which closely resemble the 
BEC wavefunction shifted in momentum space, even 
though such structures were not present in the initial 
noise which triggered the instability. This lends addi- 
tional validity to our choice of basis on which to analyze 
the atomic field dynamics. 

In conclusion, we analyze directional superradiant 



emission of atom pairs from elongated BECs by con- 
structing a simplified coupled-mode model which de- 
scribes both atom-photon and atom-atom pair genera- 
tion. Our model indicates that while still experimentally 
realizable, conditions for directional atom-atom pairs are 
significantly more restrictive than for directional atom- 
photon pairs. 
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ical Atomic and Molecular Physics at Harvard University 
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FIG. 2. Numerical simulations of the atomic density in mo- 
mentum space. Atom number is 10^, a = 10 and 5 = 10^. 
Control parameter values are (a) y = 2, (b) y = 10, (c) y = 20, 
and (d) y — 100. The single-mode diffraction angle 9d is indi- 
cated in (a), while in (b)-(d) the contour corresponds to the 
half- maximum of the gain feature of Eq. 
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